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Abstract — For nonlinear discrete time systems satisfying 
a controllability condition, we present a stability condition 
for model predictive control without stabilizing terminal con- 
straints or costs. The condition is given in terms of an analytical 
formula which can be employed in order to determine a 
prediction horizon length for which asymptotic stability or 
a performance guarantee is ensured. Based on this formula 
a sensitivity analysis with respect to the prediction and the 
possibly time varying control horizon is carried out. 

I. INTRODUCTION 

By now, model predictive control (MPC) has become a 
well-established method for optimal control of linear and 
nonlinear systems, see, e.g., [5] and [3], [20]. The method 
computes an approximate closed-loop solution to an infinite 
horizon optimal control problem in the following way: in 
each sampling interval, based on a measurement of the 
current state, a finite horizon optimal control problem is 
solved and the first element (or sometimes also more) of 
the resulting optimal control sequence is used as input for 
the next sampling interval(s). This procedure is then repeated 
iteratively. 

Due to the truncation of the infinite optimization horizon 
feasibility, stability, and suboptimality issues arise. Subopti- 
mality is naturally discussed with respect to the original infi- 
nite horizon optimal control problem, cf. [12], [17], [21], but 
there are different approaches regarding the stability issue. 
While stability can be guaranteed by introducing terminal 
point constraints [14] and [1] or Lyapunov type terminal 
costs and terminal regions [6], [16], we focus on a particular 
stability condition based on a suboptimality index introduced 
in [8], for unconsti-ained MPC — that is MPC without 
modifications such as terminal constraints and costs. Here, 
we present a closed formula for this suboptimality index. 
This enables us to carry out a detailed sensitivity analysis of 
this stability criterion with respect to the optimization and 
the control horizon, i.e., the number of elements of the finite 
horizon optimal control sequence applied at the plant. 

Typically, the length of the optimization horizon pre- 
dominantly determines the computational effort in each 
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MPC iteration and is therefore considered to be the most 
important parameter within the MPC method. However, 
suitably choosing the control horizon may lead to enhanced 
performance estimates and, thus, to significantly shorter 
optimization horizons. In particular, we prove linear growth 
of the prediction horizon for appropriately chosen control 
horizon with respect to a bound on the optimal value function 
— an estimate which improves its counterparts given in 
[7] and [22]. Furthermore, we show that MPC is ideally 
suited in order to deal with networked control systems. To 
this end, the stability proof from [8] is extended to time 
varying control horizons which allows to compensate packet 
dropouts or non-negligible delays. Here, we show that the 
corresponding stability condition is not more demanding than 
its counterpart for so called "classical" MPC for a large 
class of systems. In addition, the results in this paper lay the 
theoretical foundations for MPC algorithms safeguarded by 
performance estimates obtained for longer control horizons 
as developed in [18]. 

The paper is organized as follows. In Section III] the 
problem formulation and the required concepts of multistep 



feedback laws are given. Then, in Section III a stability con- 
dition is derived and analysed with respect to the prediction 



horizon. In the ensuing Section IV a stability theorem allow- 
ing for time varying control horizon is presented. In order 
to illustrate our results an example of a nonlinear inverted 
pendulum on a cart is considered and some conclusions are 
drawn. 

II. Problem Formulation 

In this work we consider nonlinear control systems driven 
by the dynamics 



x(n + 1) = f{x{n), u(n)) 



(1) 



where x denotes the state of the system and u the externally 
applied control. Both state and control variables are elements 
of metric spaces {X,dx) and {U,du) which represent the 
state space and the set of control values, respectively. Hence, 
our results are also applicable to discrete time dynamics 
induced by a sampled finite or infinite dimensional system. 
Additionally, state and control are subject to constraints 
which result in subsets X C X and U C [/. Given an 
initial state xq E X and a control sequence u = (u(n))„g/, 
/ = {0, 1, . . . , iV - 1} with iV e N or / = No, we denote 
the corresponding state trajectory by a;„(-) = Xu{-;xq). Due 
to the imposed constraints not all control sequences u lead 
to admissible solutions. Here, U^ [xq) denotes the set of all 
admissible control sequences u = (u(0), u(l), . . . , u{N—l)) 



of length N satisfying the conditions f{xu{n)^u{n)) e X 
and uiji) e U for n = 0, . . . , iV - 1. 

We want to stabiHze ([TJ at a controlled equilibrium x* 
and by u* we denote a control value with f{x*,u*) — 
X*. For given continuous stage costs £ : X x U -^ R^ 
satisfying i{x*,u*) ~ and £{x,u) > for all u E U for 
each x ^ x*, our goal is to find a static state feedback 
law X ^^ U which minimizes the infinite horizon cost 
Joo{x,u) = J2'^=Q^{^u{n),u{n)). Since this task is, in 
general, computationally intractable, we use model predictive 
control (MPC) instead. Within MPC the cost functional 



N 



Jn{x,u) ;= y^J{xu{n;x),u{n)) 



(2) 



is considered where A^ £ N>2 denotes the length of the pre- 
diction horizon, i.e. the prediction horizon is truncated and, 
thus, finite. The resulting control sequence itself is also finite. 
Yet, implementing parts of this sequence, shifting the predic- 
tion horizon forward in time, and iterating this procedure ad 
infinitum yields an implicitly defined control sequence on the 
infinite horizon. While typically only the first control element 
of the computed control is applied, cf. [20], the more general 
case of multistep feedback laws is considered here. Hence, 
instead of implementing only the first element at the plant 
(to = 1), m £ {1, 2, . . . , A^ — 1} elements of the computed 
control sequence u = (u(0), w(l), . . . , u{N—l)) are applied. 
As a result, the system stays in open-loop for m steps. The 
parameter m is called control horizon. 

Definition 1 (Multistep feedback law): Let A^ G N>2 and 
m £ {1, 2, . . . , A^ — 1} be given. A multistep feedback law is 
a map HN,m : X x {0, 1, . . . , m — 1} — > C/ which is applied 
according to the rule a;^„ ,„ (0; x) = x, 

Xf,{n + 1; x) = f{Xf,{n; x), fi{Xf,{<f{n); x),n- ip{n))) 

with fi — iiN,m and ip{n) := niax{fc77i|fc £ Nq, km < n}. 

For simplicity of exposition, we assume that a minimizer 
M* of (|2]i exists for each x £ X and N £ N. Particularly, 
this includes the assumption that a feasible solution exists 
for each x £ X. For methods on avoiding this feasibility 
assumption we refer to [19] or [10]. Using the existence of 
a minimizer u* £ U^ [x], we obtain the following equality 
for the optimal value function defined on a finite horizon 

Vn{x) := inf Jn{x,u) = Jn{x,u*). (3) 

Then, the MPC multistep feedback iiM.mi'i ) is defined by 
fJ-N,m{x, n) — u*{n) = u*{n; x) for n = 0, 1, . . . , ttt. — 1. In 
order to compute a performance or suboptimality index of 
the MPC feedback /i ~ fJ.N,m, we denote the costs arising 
from this feedback by 

oo 

J^{x) ■=^^ {xp^in; x),fi{Xf,{ip{n);x), n - if{n))) . 

n=0 

Notation: throughout this paper, we caU a continuous 
function p : K>o —J- M>o a class /Cco-function if it satisfies 
p(0) = 0, is strictly increasing and unbounded. 



III. Stability Condition 

In this section we derive a stability condition for MPC 
schemes without stabilizing terminal constraints or costs. To 
be more precise, we propose a sufficient condition for the 
relaxed Lyapunov inequality 

m—l 

Vat (a;^ (to; a;)) < Vn{x) - a ^ l{x^{n;x),p.{x,n)), (4) 

X £\, with a £ (0, 1] which, in turn, implies a performance 
estimate on the MPC closed-loop, cf. [12], [15]. We point out 
that the key assumption needed in this stability condition is 
always satisfied for a sufficiently large prediction horizon if 
we suppose that the optimal value function Voo (•) is bounded, 
cf. [2], [7], [13]. In particular, the formula to be deduced 
allows to easily compute, e.g., a prediction horizon for which 
stability or a desired performance estimate is guaranteed. 

Theorem 2: Let a prediction horizon A^ £ N>2 and a 
control horizon m £ {1,2,...,A^— l}be given. In addition, 
let a monotone real sequence F — (70,71572; ■•■ ,7Ar)j 
70 = 1, exist such that the inequality 

yi{x) < liViix) — 7i min l(x,u) \f x £% (5) 

ueV:f(x,u)eii 

holds for all i £ {1,2,..., A^}. Furthermore, assume that the 
suboptimality index a — aN,m given by 



Q = 1- 



N 

n (7. 

i—m-{-l 



N 

1) n (7. 

i=7V-Tra+l 



1) 



N N 

_i—m-\-l 2— m+1 



' [ N N 

1) Ul^ - n(7.-i) 

li=N~m+l i=N-m+l 

(6) 
satisfies a > 0. Then, the relaxed Lyapunov Inequality (|4]) 
holds for each x e X for the feedback law /x = HN^m and the 
corresponding MPC closed-loop satisfies the performance 

estimate 

Jlif^-{x)<-Vo.{x). (7) 



If, in addition, /Coo -functions V^ V exist such that 

ri{dx{x,x*j)<Vi{x) and Vn{x) < fj{dx{x,x*)) (S) 
hold, then the MPC closed-loop asymptotically converges to 

X*. 

Proof: We sketch the main ideas of the proof and refer 
for details to [11] for the main part and to [24] for the 
adaptation to our more general setting. 

Using Bellman's principle of optimality and Condition Q 
in order to derive conditions on an open-loop optimal trajec- 
tory allows to propose the following optimization problem 
whose solution yields a guaranteed degree of suboptimality 
a for the relaxed Lyapunov Inequality Q: 



inf 

Xq,...,\j^^I,1' 



Z^n=0 ^n 



E 






subject to the constraints 

JV-l 

^ A„ < 7JV-fe-Afc, k^0,...,N-2, 



n^k 



j-i 



- 22 -^"^ 



< 



In-] ■ Aj+„i, 



J=0, 



,iV-m-l, 



n=0 



and Aq, . . . , \n-i, i^ > 0- Here, we used the abbreviations 
A„ := £{xu*{n),u*{n)) for a minimizer w* € lA^ {x) of (|2]) 
and z/ := Viq{xu*{m)). 

Within this problem, the constraints represent estimates 
obtained by using (|5]) directly or first following an optimal 
trajectory and, then, making use of (|5|. In the next step, this 
optimization problem is reformulated as a linear program. 
Then, neglecting some of the imposed inequalities leads to a 
relaxed linear program whose solution is given by Formula 
(|6]l. Hence, a from Formula (|6]) is a lower bound for the 
relaxed Lyapunov Inequality Q. 

If the submultiplicativity condition 



A„r • A„,r > A„+,„r with /\iV:=-ii 



(9) 



is satisfied for all n.ra ^1^ with n + m < N and the given 
sequence F, Formula (|6]l actually solves the non-relaxed 
problem and, thus, characterizes the desired performance 
bound even better. Otherwise, solving the non-relaxed prob- 
lem may further improve the suboptimality bound a. ■ 
Remark 3: The main assumption in Theorem l2] is In- 
equality (|5]l which is also used in [7], [22]. However, the 
performance estimates deduced in these references are more 
conservative in comparison to the presented technique, cf. 
[23]. The controllability condition used in [11], i.e. existence 
of a sequence (c„)„gNo C M>o such that for each state x E X 
an open-loop control u^ E U°°{x) exists satisfying 



l{xu^{n;x),Ux{n)) < CnVi{x) 



(10) 



implies Inequality (|5]l with -fi :— X^n^o "^" ^^'- l^ads, in 
general, to more conservative estimates, cf. [23]. Note that 
the methodology proposed in [8] allows to use sequences 
r = (7i)ieN depending on the state. Furthermore, we 
emphasize that a suitable choice of the stage costs may 
lead to smaller constants 7^, i G {1,2,...,A^} and, thus, 
to improved guaranteed performance, cf. [4] for an example 
dealing with a semilinear parabolic PDF. 

We like to mention that Theorem |2] can be extended to 
the setting in which an additional weight on the final term 
is incorporated in the MPC cost functional, i.e. 



Jn{xo,u) 



N~2 

E 

n=0 



l{xu{n),u(n))+ujl{xu{N-l),u(N~l)) 



with i:j > 1, cf. [11, Section 5]. 

The availability of an explicit formula facilitates the analy- 
sis of the performance estimate aN.m and, thus, allows to 
draw some conclusions. The first one, stated formally in the 
Corollary HI below, is that MPC without stabilizing terminal 
constraints or costs approximates the optimal achievable 
performance on the infinite horizon arbitrarily well for a 



sufficiently large prediction horizon N — independently of 
the chosen control horizon m. For the proof, the concept of 
an equivalent sequence given in [23] is employed. Then, the 
argumentation presented in [11, Corollary 6.1] can be used 
in order to conclude the assertion. 

Corollary 4: Let the controllability Condition (|5]l be sat- 
isfied for a monotone bounded sequence F = (7i)ieN- 
Furthermore, let a control horizon ttt, e N be given. Then, 
the suboptimality estimate ajv.m. ^ > max{2,TO + l}, from 
Formula (|6]) converges to one for N approaching infinity, i.e. 
limAr_j.oo aN,m = 1- If, in addition. Condition (|8]l holds, the 
MPC closed-loop is asymptotically stable. 

In order to further elaborate the benefit of Formula (|6|, 
the following example is considered. 

Example 5: Let an exponentially decaying function 
/3(r, n) = C(j"r be given. Then, for each prediction horizon 
N £ {2, 4, 8, 16}, we determine all parameter combinations 
{C,cr) e M>i X (0,1) such that the stabiHty condition 
ajv,! > holds with 7, = CX^^lgCr", cf. Fig. Jll Note 
that this setting corresponds to assuming Condition {T0[ with 
c„ = Ca". 
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Fig. 1. Parameter pairs (C, o") for function /3(r, n) = Ccr"r such tliat tlie 
corresponding performance bound satisfies the stability condition ajv 1 > 
depending on the prediction horizon A'^. 



We point out that Fig. [T] shows the different influence of 
the overshoot C and the decay rate a. Indeed, the figure 
indicates that for given N > 2 and a £ (0, 1) stability always 
holds if the overshoot C > 1 is sufficiently small. However, 
for given N > 2 and overshoot C > the stability condition 
may be violated regardless of how a £ (0, 1) is chosen. This 
observation can be proved rigorously using Formula (|6|, cf. 
[11, Proposition 6.2]. 

Secondly, Theorem l2] allows to deduce asymptotic esti- 
mates on the minimal prediction horizon length N for which 
the stability condition ajv.m > 0, m £ {1, 2, . . . , A^ — 1}, 
holds — depending on the sequence F ~ {li)ien from 
Condition (|5]l. Here, one has to keep in mind that the 
prediction horizon N predominantly determines the required 
computation time in order to solve the finite horizon opti- 
mization problem in each iteration of an MPC algorithm. 



The next proposition uses a special version of Inequality 
(|5]l in which the 7^ are independent of i. It can be checked, 
for instance, using an upper bound for the optimal value 
function Voo, cf. [11, Section 6] for a proof. 

Proposition 6: Let Condition (|5]l be satisfied with F = 
(7i)ieN with 7i = AI for all i e NP] Then, asymptotic 
stability of the MPC closed-loop is guaranteed if, 

• for m ~ I, the following condition on the optimization 
horizon is satisfied 

N>2 + ^"("^^ - ^) (11) 

- +ln(A/)-ln(A/-l) ^ ' 

and, thus, the minimal stabilizing prediction horizon 

N := min{7V : N e N>2 and ajv,™ > 0}. (12) 

grows asymptotically like Mln(M) as M — >■ 00, 

• for 771 = [N/2\, one of the following inequalities holds 



N > 



N > 



21n(2) 



ln(M) - lii(Af - 1) ' 

'2M-1\ 1 /2M-1N 



N even (13) 



ln(^ 



M 



-)H 



Af-l 



A^odd. 



(14) 



\n{M) - ln(A/ - 1) ' 

In this case, the minimal stabilizing Horizon ( [T2| ) grows 
asymptotically like 21n(2)A/ as M — ?► 00. 
By a monotonicity argument, the estimates from this propo- 
sition also apply to each sequence T — (7i)igN which is 
bounded by M. 




Fig. 2. Minimal prediction liorizon A'^ for wliicli stability is guaranteed by 
Theorem 12] supposing controllability Condition J5) with (7i)igN. 7i = ^'1 
for all i eN. 

The conclusions of Proposition l6] are twofold: First, nu- 
merical observations from [8] are confirmed and the cor- 
responding parameters are precisely determined. Secondly, 
we emphasize the linear growth of the minimal stabilizing 
prediction horizon for m = [N/2\ . Hence, the growth for 
larger control horizons is much slower than for MPC with 
control horizon m = 1, cf. Fig.l2] This fact will be exploited 

'Note that the value of 71 is not taken into account in the computation 
of Q^]V,m from Formula J6}. Indeed, 72 is the first value contributing to the 
corresponding suboptimality index. 



in the following section for both networked control systems 
and for "classical" MPC by designing suitable algorithms. 

IV. Time Varying Control Horizon 

In the previous section a stability condition was derived 
which can also be used in order to ensure a guaranteed 
performance of the MPC closed-loop. As Proposition |6] 
already indicated, employing larger control horizons may 
improve the corresponding estimates on the required predic- 
tion horizon length for which stability can be guaranteed. 
The following proposition states further properties of the 
suboptimality Bounds (|6|. 

Proposition 7: Suppose that Condition (|5]l holds with T = 
(7^)^eN, 7^ -^ C Y^'}^ <j'\ Here, C > 1 and a e (0,1) 
denote overshoot and decay rate of a system which is 
exponentially controllable in terms of the stage costs. Then, 
the performance Estimate (|7]i has the properties: 

• symmetry, that is aN.m — oiN.N-m, and 

• monotonicity, i.e. a^q m+i > o-n m for all m G 
{l,2,...,LiV/2j-l}.^ 

As a consequence, a^.m > ctN,i holds for all m G 
{1,2,...,A^— 1} and, in particular, the stability condition 
ctN,m > holds for arbitrary control horizon m > 2 if it is 
satisfied for ni = 1. 

Proof: Symmetry follows directly from Formula (|6|. 
Contrary to this, showing the claimed monotonicity prop- 
erties requires a more elaborate technique, cf. [11, Section 
7] for a detailed proof. ■ 

Proposition IT] can be exploited in various ways. For 
instance, in networked control systems the fact ajv.m > cun.i 
for all m £ {l,2,...,iV — 1} can be used in order to 
conclude stability of a compensation based networked MPC 
scheme in the presence of packet dropouts or non-negligible 
delays. The compensation strategy is straightforward: Instead 
of sending only one control element across the network, an 
entire sequence is transmitted and buffered at the actuator If 
a packet is lost or arrives too late — that is the packet has 
not been received by the actuator by the time the first control 
element of this sequence has to be implemented — the 
succeeding element of the current sequence is implemented 
at the plant which corresponds to incrementing the control 
horizon m. Since it is a priori unknown when and if the next 
package and, thus, the next sequence of control values arrives 
at the actuator, the control horizon has to be time varying. 
Using Theorem [8] stability can nevertheless be concluded. 

In order to formulate this assertion in a mathematically 
precise way, the following notation is needed: Let to* £ 
{2, . . . , iV — 1} be an upper bound for the maximal number 
of elements of the computed control sequence to be imple- 
mented. Then, the transmission times are given by a sequence 
of control horizons M = {mk)i^^j^^ with to* > to.^ > 1. 
Consequently, in between the fcth and the (fc + l)st update 
of the contol law the system stays in open-loop for nik 
steps. Here, we denote the update time instants by cr(fc) := 
Si=o "^i while (p{n) := max{cr(A:) | k G Nq, a{k) < n} 
maps the time instant n € No to the last update time 
instant. The corresponding control law is denoted by hn,m- 



Illustrating these new elements, a control sequence is a 
sequence 

^i{Xf_,{a{k);x),0),... ,^i{x^{a{k)\x),mk - 1), 
^{x^{a{k + l);x),Q),... 

with /i = y.NM- 

Theorem 8: Suppose that a multistep feedback law 
AiAr,,„* : X X {0, . . . ,m* - 1} ^ f/, m* < iV - 1, and 
a function Vn : X -^ R^ are given. If, for each control 
horizon m e {1, 2, . . . , m*} and each x £ X, we have 



-1 

E 

n=0 



V"(a;) - V"(a:^(TO; x))>a} l{Xf,{n\ x),fi{x, n)) (15) 



with /i = fJ.N.m* for some a G (0, 1], then the estimate 
avj{x) < aV^"'" {x) < VNix) holds for all a; e X and all 
M = (mfc)fcgNo satisfying irik < m*, k G Nq. If, in addition. 
Condition ([HJ is satisfied for VAr(), asymptotic stabiUty of 
the MPC closed-loop is ensured. 

Theorem [8] generalizes its counterpart [8, Theorem 5.2] to 
time varying control horizon. To this end, the value function 
VAr(-) was used as a common Lyapunov function, cf. [11, 
Theorem 4.2]. In order to verify the required assumptions of 
Theorem [8] our stability condition has to hold for different 
control horizons m, i.e. for each m £ {1,2,..., m*} which 
can be checked by Theorem |2] However, e.g. for an ex- 
ponentially controllable system. Proposition [7] automatically 
ensures this condition if it is satisfied for m — 1. Hence, 
the stability condition for time varying control horizons 
remains the same as for MPC with m = I. Furthermore, 
we like to point out that increasing the control horizon often 
enhances the proposed suboptimality bound significantly. In 
particular, this improvement may lead to a stability guarantee 
by aN,rn > although this conclusion cannot be drawn for 
m ~ 1 (ajy,! < 0), cf. Fig. [S] and the numerical results 
shown in Section |V] 




Fig. 3. Let Condition J5j be satisfied for the sequence F = (7i)igN with 
7i = C^l^Q it", C = 3, ct = 2/3. Then, the smallest hoiizon for which 
o^ AT, 771 ^ and, thus, stability is guaranteed by Theoreml2|is N = 12. For 
m = 1, even a prediction horizon of length A'^ = 18 is required. 

Another way to use Proposition l7] is described in [18]. 
There, an algorithm is constructed which employs larger 



control horizons in order to guarantee a desired performance 
bound. Then, based on an evaluation of the relaxed Lyapunov 
Inequality ( [T5| ), the MPC loop is closed as often as possible 
performing a new MPC optimization. This procedure often 
leads to MPC with m — 1, however, safeguarded by the 
fact that the desired performance can always be ensured — 
if necessary — by enlarging m, cf. Fig. [3] The observed 
improvement be explained as follows: Checking the relaxed 
Lyapunov Inequality ([TSj at each time instant is a sufficient 
but not a necessary condition for ^T5\ to hold for m > 1, i.e. 
larger control horizons lead to less restrictive conditions. 

V. Example 

We illustrate our results by computing the a-values from 
the relaxed Lyapunov Inequality ([TSj along simulated tra- 
jectories in order to compare them with our theoretical 
findings. We consider the sampled-data implementation of 
the nonlinear inverted pendulum on a cart given by the 
dynamics 

Xi{t) = X2{t) 

X2{t) = -j Sin{xi{t) + ^) - -yX2{t)\x2{t)\ 



— u{t) cos{xi{t) + it) 

Xsit) = X4{t) 
X4{t) = u{t) 



kBSgn{x2{t)) 



where g = 9.81, I = 10 and kji = k^ = 0.01 denote 
the gravitation constant, the length of the pendulum and 
the air as well as the rotational friction terms, respectively. 
Hence, the discrete time dynamics is defined by x{n ~\- 
1) = $(T; a:(n),M(n)). Here, $(T; a;(n),u(n)) represents 
the solution of the considered differential equation emanating 
from x{n) with constant control u{t) = u{n), t e [0,T) at 
time T. The goal of our control strategy is to stabilize the 
upright position x* — (0, 0, 0, 0). To this end, we impose the 
stage cost 



i{x{n),u{n)) -.^ / l{<^{t;x{n),u{n)),u{t))dt 
Jo 

with l(x,u) given by 

lO-^u^ + (3.51 sin^ xi + 4.82x2 sinxi + 2.31a;^ 

+ 0.01 xl + 2 ((1 - cosa;i)(l + cos^ X2)f + O.lxl) 

with sampling time T = 0.05 and prediction horizon N = 
70. Within our computations, we set the tolerance level 
of the optimization routine and the error tolerance of the 
differential equation solver to 10^^ and 10^^, respectively. 
Due to the 2tt periodicity of the stage cost £, we limited the 
state component xi to the interval [— 27r + 0.01, 27r — 0.01] in 
order to exclude all equilibria of ^ different from x*. All other 
state components as well as the control are unconstrained. 
For our simulations, we used the grid of initial values 



g ■^{xeR^lBiei-i.O,!}^ 



0.05i} 



with X = (tt + 1.4, 0, 0, 0)^ and computed the suboptimality 
degree ayo.m for constant control horizons nik = rn along 
the MPC closed loop. 

Here, we used a startup sequence of 20 MPC steps with 
TO = 1 to compensate for numerical problems within the 
underlying SQP method. The startup allowed us to compute 
an initial guess of the optimal open-loop control close to the 
optimum. During our simulations, we were able to achieve 
practical stability only, a fact we compensated within our 
calculations by introducing a truncation region of the stage 
cost £ using the constant e — 10^"'. The idea of this cut 
is to take both practical stability regions, that is small areas 
around the target in which no convergence can be expected, 
and numerical errors into account, cf. [9, Theorem 21] for 
details. The values of aN,m are computed along the closed- 
loop trajectory via 

Oat,™ = min inf aN,m{n;xo) (16) 

xqEG nS{ri|3/cGNo:n=fcm} 

with local degree of suboptimality aAr,m(") given by 

V/v(a;p(n; xq)) - VAr(.T^(n + to; xq)) 



aN,min;xo) = 



rn— 1 

J2 Wxfiin + k; xo),ii{k; a;^(n; Xq))) - e) 

fc=0 



with /i = ^J'N,m if the denominator of the right hand side 
is strictly positive and aN,„i{n;xo) = 1 otherwise. Note 
that ttTv.m may still become negative if the value function 
increases along the closed-loop. 

In Fig. HI ajv.m-values according to ( fTSj l are shown for a 
variety of control horizons to using the optimization horizon 
iV = 70. 
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Fig. 4. Approximation of «7o.m, »" £ {1, . . . , 70} for the nonlinear 
inverted pendulum. 



While for to < 11 stability of the closed loop cannot 
be guaranteed, we obtain avo.m > for to e [12,47]. For 
r?i > 48 the values of ayo .„ are decaying rapidly which may 
be the result of numerical problems. 



VI. Conclusions 

We presented a stability condition for MPC without termi- 
nal constraints or Lyapunov type terminal costs for nonlinear 
discrete time systems, which allows to determine a prediction 
horizon length for which asymptotic stability or a desired 
guaranteed performance is ensured. Furthermore, we investi- 
gated the influence of the prediction and the control horizon 
on this condition. Suitably choosing the control horizon 
leads to linear growth of the prediction horizon in terms 
of the assumed controllability condition. As a consequence, 
since the prediction horizon predominantly determines the 
computational costs, computing times can be reduced. In 
addition, a stability theorem for time varying control horizons 
was derived. Using symmetry and monotonicity properties, 
we showed that no additional assumptions were needed in 
comparison to "classical" MPC. 
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